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ABSTRACT 

We announce the discovery of a low-mass planet orbiting the super metal-rich KOV star 
HD77338 as part of our on-going Calan-Hertfordshire Extrasolar Planet Search. The best fit 
planet solution has an orbital period of 5.7361±0.0015 days and with a radial velocity semi- 
amplitude of only 5.96±1.74 ms" 1 , we find a minimum mass of 15.9ig;3 M e- The best fit eccen- 
tricity from this solution is 0.091q'q9, and we find agreement for this data set using a Bayesian 
analysis and a periodogram an alysis. We measure a metallicity for the star of +0.35±0.06 dex, 



whereas another recent work ( Trevisan et al.l 120111 ) finds -f0.47±0.05 dex. Thus HD773386 is 



one of the most metal-rich planet host stars known and the most metal-rich star hosting a sub- 
Neptune mass planet. We searched for a transit signature of HD773386 but none was detected. 
We also highlight an emerging trend where metallicity and mass seem to correlate at very low 
masses, a discovery that would be in agreement with the core accretion model of planet forma- 
tion. The trend appears to show that for Neptune-mass planets and below, higher masses are 
preferred when the host star is more metal-rich. Also a lower boundary is apparent in the super 
metal-rich regime where there are no very low-mass planets yet discovered in comparison to the 
sub-solar metallicity regime. A Monte Carlo analysis shows that this low-mass planet desert is 
statistically significant with the current sample of 36 planets at ~4.5tr level. In addition, results 
from Kepler strengthen the claim for this paucity of the lowest-mass planets in super metal-rich 
systems. Finally, this discovery adds to the growing population of low-mass planets around low- 
mass and metal-rich stars and shows that very low-mass planets can now be discovered with a 
relatively small number of data points using stable instrumentation. 

Subject headings: stars: fundamental parameters — stars: (HD77338) — stars: rotation — (stars:) 
planetary systems 
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1. Introduction 

One of the first correlations announced between 
two parameters that dealt with exoplanets and 
their host stars was the abunda nce of heavy ele- 
ments in the stellar atmospheres. iGonzalea (J1997) 
first noted that all three exoplanet host stars 
known at that time were over abundant in iron. 
The paper indicated that the metallicity of these 
stars were all in the super solar regime. This fea- 
ture has since been stu died by a number of au - 
thors, most notably by iFischer fc Valentil (j2005l ) 
who defined a relationship between the host star 
metallicity ([Fe/H]) and the probability of a star 
hosting a gas giant planet. 

This metallicity bias was one of the major fea- 
tures that helped to confirm that the core accre- 
tion scenario of planet forma tion (e.g. Ilda fc Lin 



2004at iMordasim et al.l 120091), coupled w ith plan 



etary migration (|Lin &: Papaloizoulll986l) , appears 
to be the dominant mechanism to build these plan- 
etary systems. A further validation of the core ac- 
cretion mechanism again comes from stellar metal- 
licity since it appears that low-mass Neptunes and 
super-Earths are not predominant ly found around 



2006). 



metal-rich stars (|Udrv et al 

The Calan-Hcrtfordshire Extrasolar Planet 
Search (CHEPS) is a program to mainly moni- 
tor metal-rich stars to primarily hunt for short 
period planets across a wide range of masses and 
better constrain the statistics of planets around 
these stars, whilst also searching for planetary 
transits of brig ht and nearby stars in the southern 
hemisphere. In I Jenkins et al.l (|2008l ) we discussed 
the sample selection for the CHEPS, using pre- 
vious observations made with the ESO-FEROS 
instrument to measure the chromospheric activ- 
ity and metallicity of a sample of a few hundred 
stars. This sample ha s recently been incr eased 
using data discussed in Ijenkins et al.l ( 2011 ) and 
current metallicity analyses of these stars is still 
ongoing using our new methods for measuring ac- 
curate atomic abundances in stars like the Sun 
(|Pavlenko et al.ll2012l : Ijenkins et al.ll2012t ). 

In Ijenkins et al.1 ((20091) we published the dis- 



1 Email: jjenkins@das.uchile.cl. Based on observations 
collected at the La Silla Paranal Observatory, ESO (Chile) 
with the HARPS spectrograph on the ESO 3.6m telescope, 
under the program IDs 079.C-0927, 081.C-0148, 087.C- 
0368, and 088.C-0662. 



covery of an eccentric brown dwarf, or ex- 
treme Jovian planet, orbiting the metal-rich star 
HD191760, along with new orbits for three other 
recently discovered southern hemisphere metal- 
rich planets. The CHEPS data has mostly 
been obtained using the HARPS instrument 
but significant observing time on CORALIE has 
also recently been acquired to pre-select inter- 
esting new targets fo r HARPS follow-up (see 
Jenkins fc Jordan 2010l ). 



2. Observations and Reduction 

All the radial velocity data we present in this 
work were obtained usi ng the HARPS instru- 
ment (jMavor et al.l l2003h . From empirical data, 
HARPS has been shown to be r adial velocity sta - 
ble down to the sub- lms -1 level ( Pepe et al1l2011 ). 
HARPS itself is a fiber-fed cross-dispersed cchcllc 
spectrograph that employs two 1" fibers, which in 
simultaneous Thorium- Argon (ThAr) mode places 
one fibre on the source (HD77338 in this case) and 
one feeds the calibration ThAr lamp to monitor 
the instrumental drift. The calibration lamp is 
useful for sub-lms _1 nightly stability as the instru- 
ment itself is stable at the lms _1 level throughout 
the course of a single night. 

We discuss 32 velocity measurements from 
HARPS in this current work for the star HD77338 
that are shown in Table [1] We follow the proce- 
dure of observing our targets over multipl e nights 
in any single run (see Jenkins et al. 20091 ). which 
are usually 3 to 5 days in duration, and with a to- 
tal integration time of 15 minutes to average over 
the st rongest p-mode oscil lations in these types of 



stars (jO'Toole et al.ll2009h 



Our latest data were reduced and velocities ex- 
tracted in near real time at the telescope using 
the most up to d ate version of the HARPS-DRS 
(jPepe et al.ll2004l ). However, our older data, which 
were acquired before the latest version of the soft- 
ware was installed on site, were re-reduced offline 
in the ESO Vitacura offices using the latest version 
of the DRS. This software performs all functions 
necessary to fully reduce and analyse a radial ve- 
locity timeseries, returning barycentric corrected 
velocities, uncertainties, and a suite of diagnostics 
to test the line stability of the data set. All ac- 
tivities we analyse were extracted using a custom 
code developed following our previous experience 
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with this type of chromospheric act ivity analysis 
(e.g. IJenkins et aDl2006l 120081 l201ll ). 



HD77338 



The star HD77338 is classified as KOIV in 



Table 1: HARPS radial velocities for HD77338. 



BJD 



RV 



Pint a tot 



2453358 


7928045 


-18.11 


1 


66 


2 


67 


2453359 


7935070 


-7.53 


1 


89 


2 


83 


2453360 


8256229 


-4.32 


1 


61 


2 


65 


2454058 


8168860 


-15.45 





49 


2 


16 


2454098 


6739885 


-19.44 





71 


2 


22 


2454103 


7481653 


-16.52 





73 


2 


15 


2454124 


6791307 


-9.52 





46 


2 


15 


2454125 


7692682 


-15.18 





55 


2 


17 


2454133 


7288993 


-12.63 





55 


2 


17 


2454161 


6324334 


-15.91 


1 


09 


2 


37 


2454161 


6593658 


-18.24 





46 


2 


15 


2454162 


6659105 


-12.47 





49 


2 


16 


2454365 


8855969 


-10.19 





74 


2 


23 


2454579 


6605011 


-22.23 





78 


2 


24 


2454580 


6290112 


-17.33 





72 


2 


22 


2454581 


5835519 


-7.64 





65 


2 


20 


2454725 


8665876 


-10.20 





82 


2 


25 


2454726 


8663334 


-10.79 





72 


2 


22 


2454726 


9022721 


-8.34 





58 


2 


18 


2454727 


8587883 


-14.39 


2 


09 


2 


96 


2455650 


5916958 


-10.35 





65 


2 


20 


2455651 


6037866 


-13.69 





75 


2 


23 


2455883 


7133405 


-11.57 





76 


2 


24 


2455883 


8696951 


-12.36 





59 


2 


18 


2455885 


7293282 


-12.72 


1 


03 


2 


34 


2455885 


8583195 


-10.69 


1 


02 


2 


33 


2455992 


5489188 


-9.95 





48 


2 


16 


2455992 


7888091 


-7.77 





65 


2 


20 


2455993 


5357474 


-5.27 





48 


2 


15 


2455993 


7774660 


-6.55 





56 


2 


17 


2455994 


5078761 


-7.82 





60 


2 


18 


2455994 


7518050 


-10.03 





68 


2 


21 



The ai n t and atot are the uncertainties on the 
points before and after adding the estimated 
activity jitter value shown in Table [3] in 
quadrature. 



the Hipparcos main catalogue (jPerrvman et al 
1997T1 and has a parallax of 24.54±1.06 mas 



(jvan Leeuwen 2007 ). giving rise to a distance of 
40.75±1.76 pc. It has an apparent U-band magni- 
tude of 8.63, meaning its distance modulus is 3.05, 
and the Johnson B — V colour index is 0.833. 

The photometric and astrometric analysis 
means we can place HD77338 on a HR-diagram 
and compare its position to that of model isochrones 
and isomass tracks. However, in order to do this 
properly, we must first have a handle on the metal- 
licity. Given that the star is fairly nearby and 
bright, there have been a number of attempts to 
pin down its metallicity. 

Using FEROS spectroscopy we had previ- 
ously measured the [Fe/H] of HD 77338 to be 
+0.25±0.07 dex (| Jenkins et al.ll2008l ) clearly clar- 
ifying that this star is indeed super metal-rich. 

Since we now have timeseries data from HARPS 
we have also measured the metallicity, and other 
atomic abundances using this spectra. We nor- 
malised each of the ID spectra and combined 
all the residual fluxes to gain a high S/N and 
high resolution single HD77338 optical spectrum. 
We then applied our new method for measuring 
self-consistent abundances from many optical ab- 
sorption lines for Fe 1 and Fe 11 independently 
(|Pavlenko et al.ll2012h and find a logN(Fe i) of - 
4.08±0.02 dex and a logN(Fe 11) of -4.07±0.06 dex. 
Although in agreement, the Fe 11 abundance has 
a larger uncertainty due to the lower number of 
lines used to determine this value. The final abun- 
dances for iron and other elements, along with the 
solar values and the number of lines used in the 
abundance determination for each are shown in 
Tabled 

Taking these two abundances an d subtracting 
off th e solar abundance value of -4.42 (iPavlenko et al 
20121 ) we find a [Fe/H] of +0.35±0.06 dex. This 
value is in agreement with our previous value to 
within the uncertainties but is even more metal- 
rich. It is also in good agreement with othe r au- 
thors, for example Barbuv fe Grenor] ( 1990h find 
a [Fe/H] of +0.22 dex. iTavlorl (j2005l ) find a value 
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Table 2: Elemental abundances for HD77338. 



Ion 


logJV(X) 


log7V(X) 


[X/H] 


[X/Fe] 


N, 


Si I 


-4.00±0.06 


-4.45i0.05 1 


0.45 


0.10 


13 


Si II 


-3.78±0.03 


-4.45i0.05 1 


0.67 


0.32 


2 


Ca I 


-5.25±0.02 


-5.64i0.02 1 


0.39 


0.04 


12 


Ti I 


-6.67±0.03 


-6.98i0.06 1 


0.31 


-0.04 


29 


Ti II 


-6.76±0.05 


-6.98i0.06 1 


0.22 


-0.13 


27 


Cr I 


-6.04±0.03 


-6.33i0.03 1 


0.29 


-0.06 


38 


Fe I 


-4.08±0.02 


-4.42i0.03 2 


0.34 


-0.01 


130 


Fe II 


-4.07±0.06 


-4.42i0.03 2 


0.35 


0.00 


18 


a 


-5.74i0.09 


-6.08i0.11 


0.34 


-0.01 


119 



1 iGrevesse &; Sauvall (|l998l ). 2 iPavlenko et al.l (|2012h : a = mean of Si I. Ca I. Ti I, Ti II, and Cr I. 



of +0.283i0.051 dex. ICenarro et al.l (l2007h quote 



a valu e of +0.25i0.10 dex, and iTrevisan et al 
(|201ll ) find i0.47i0.05 dex. Finally, we can be 
sure from these analyses that HD77338 is a super 
metal-rich star. 

In Fig. [T]we show the position of HD77338 on 
a HR-diagram (cross hairs), with \og(T e ff) and 
T e ff on the lower and upper x-axes, respectively, 
and log(L/L ) and L/L Q on the left and right 
y-axes, respectively. Yonsie-Yale (Y2) evolu tion 



ary tracks are shown ( Demarque et al. 2004). and 



the analysis method we perform is explained in 



Jenkins et al 



(200§). The important point here 
is that the solid curves represent three masses of 
0.8, 0.9, and 1.0M Q , from bottom to top, and the 
dashed curves show the tracks at a higher reso- 
lution of 0.01M Q within those three mass regimes. 
This analysis yields a mass of 0.93i0.05MQ, a ra- 
dius of 0.88i0.04R Q , a log g of 4.52i0.06 dex, 
and an age of 3.97tg;^ Gyrs for HD77338. We 
note that HD77338 is not a subgiant star, as la- 
beled in Hipparcos, but is a typical metal-rich K- 
dwarf star. All calculated values for HD77338 are 
shown in Table [3] 

4. Doppler Analysis 

The radial velocity timeseries of data for 
HD77338 show evidence for a signal with a very 
short period of only 5.74 days. The signal is ap- 
parent using two separate analysis methods that 
we describe here. 

4.1. Periodogram Analysis 

A periodogram search for strong and stable 
frequencies in the radial velocity data set for 



HD77338 reveal a single prominent peak. In Fig. [2] 
(upper plot) we show the periodogram for the 
HD77338 data and the strongest peak that relates 
to the 5.74 day Doppler signal is clear. At present 
there are still some small peaks around this pe- 
riod, but not with a power close to approaching 
the 5.74 day period. 

Two other peaks are seen in the periodogram, 
at periods of ~4.5 days and 1.2 days. We have 
on occasion observed this star multiple times per 
night to help break this aliasing, but it will re- 
quire further observing runs to suppress the power 
around these frequencies. In fact in our final three 
night run, we observed the star twice per night 
and this significantly reduced the power in the 
1.2 day signal, indicating this is probably a sam- 
pling frequency. We also note that by looking 
for power in the sampling, we find some power 
at periods close to one day, further validating the 
claim that the 1.2 day peak might be a sampling 
frequency and not a genuine Doppler frequency. 
The strongest peaks in the sampling frequencies 
are around 35 days, no significant power is found 
close to the 5.7 day period of our planet candidate 
signature. 

We ran a false alarm probability (FAP) test 
by scrambling the velocities with replacement, but 
retaining the time-stamps from our observations, 
to test how strong the signal we detect is, in 
comparison to its emergence being solely through 
random noise fluctuations. We ran this boot- 
strap simulation 10000 time s, similar to that in 



Anglada-Esc ude et al.l ( 2012 ). and checked how of- 



ten the peak power of the scrambled data was 
higher than or equal to the power of the observed 
data. This test revealed that 0.01% of the time 
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Table 3: Stellar parameters for HD77338. 





HD77338 




R A T9000 Hrm-c^ 

IViT. JZiUUU I 11. 111. o I 


flQ n m m 1 9 s AQA 


PArrvmnn pt ftl f 1 QQ7^ 
r lii y niciii l i cii. ^ lyj i i 


L/CL JZUUU IQ.IIl.b) 


-ZO Ol O f .^z 


jr^ciivniaii ei ai. ^iyy/ I 


Spectral Type 


ivUl V 


rcn ymdn ci ai. ^iyy/ j 


Spectral Type 


WAV 


inih vvorK 




u.ooo 


l-'ci i't-t rm -'i t i ol 11 QQ ; 1 

i i y nidii ll di. ^ryyi f 




D-UU 


Pprnrmfln fll ( 1 QQ7 1 
± Lily iiidii l l cii. ^ i y y ( ( 


distance (pc) 


/in 7£_1_1 7£ 
4U. ( o±l. 1 


van Leeuwen (2007) 


loe\n -trt^ 

i'-'S- 1 ^ rlls. 


-5.05 


Jenkins et al. (2011) 


Hipparcos iV b s 


172 


Perrvman et al. (1997) 


Hipparcos a 


0.013 


Perrvman et al. (1997) 


AM V 


0.224 


Jenkins et al. (2011) 


L,/L Q 


0.55±0.02 


This Work 


M*/M e 


0.93±0.05 


This Work 


R+/Rq 


0.88±0.04 


This Work 


Teff (K) 


5370±82 


This Work 


[Fe/H] 


0.25±0.07 


Jenkins et al. (2008) 


[Fe/H] 


0.35±0.06 


This Work 


log 3 


4.52±0.06 


This Work 


U,V,W (km/s) 


39.1,-27.7,-24.8 


Jenkins et al. (2011) 


-Prot,R'HK (days) 


45 


This Work 


-Prot,v sin i (days) 


21.7 


This Work 


Vrotfl'HK (km/s) 


1.1 


This Work 


wsim(km/s) 


2.8±1.5 


Jenkins et al. (2011) 


l>sini(km/s) 


2.33±0.05 


This Work 


Age^HK (Gyrs) 


7.5 


This Work 


Age (Gyrs) 


o 07+4-23 
°- a ' -3.00 


This Work 


Jitter - Smw ( m / s ) 


2.10 


Isaacson & Fischer (2010) 


Jitter - fit (m/s) 


1.59 


This Work 
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T EFF (K) 

4720 4960 5200 5440 5680 5920 
0.0 \ 



-0.1 - 



-0.2 - 



-0.3 



-0.4 



-0.5 




3.680 3.696 3.712 3.728 3.744 3.760 

iog(T ef( ) 

Fig. 1. — Y2 evolutionary tracks for masses of 1.0, 
0.9, and 0.8MQare shown by the solid curves from 
top to bottom. The dashed curves are the higher 
resolution analysis tracks with a step of ±0.01 Mq. 
The cross hairs mark the position of HD77338. 



* 6 




10 100 1000 10000 

Orbital Period (Days) 




10 100 1000 10000 

Orbital Period (Days) 

Fig. 2. — Lomb-Scargle periodogram of the 
HD77338 radial velocity data set is shown in the 
top plot, with the same analysis performed on the 
residuals to the 1-planet fit shown in the lower 
plot. The dashed line boundaries mark the 1% 
FAP limit. 



the peak power was larger than the observed data 
set. This gives us confidence that the signal we 
detect is robust. The 1% FAP boundary is shown 
for reference on the plots. We also note that the 
FAP registers as 0.02% if we bin the data points 
that were observed on the same night. 

We make it clear that the signal peak was al- 
ready significant with only 27 data points i.e. be- 
fore our final three night observing run. At that 
point the FAP registered as 0.02%, and this led us 
to obtain two radial velocities per night in our last 
three night run, since we should sample more than 
half of the planetary candidate orbit in that time. 
Since the signal power increased significantly from 
the addition of these six velocities, we believe it is 
a genuine Doppler reflex motion of the star and 
not related to a harmonic of our sampling strat- 
egy- 

The lower plot in Fig. [2] shows the periodogram 
of our residuals after subtraction of our best fit 
solution. There is no significant power left in the 
data after this single planet fit, an expected result 
given the relatively small number of data points we 
have acquired at present. This also highlights how 
well the Keplerian fit models the current data set 
as there is no residual power except white noise. 
The strongest frequency is found around 9.7 days, 
which could be an emerging planetary signature 
given the high frequency of multi-planet low-mass 
systems. However, much more data is required to 
confirm if this frequency is truly a Doppler induced 
signal. 

4.2. Bayesian Search 

We also run a Bayesian search for signals by as- 
suming that 0,1, 2, ...etc Keplerian signals best de- 
scribe the data we have and find the best fit that 
conforms to our assigned probability threshold. 
We analysed the radial velocities using Markov 
Chain Monte Carlo sampling of the parameter 
space of each model with a different number of sig- 
nals. These samplings were p erformed using the 



adaptive Metropol is algorithm (|Haario et al.ll2001 



iTuomi et all l201ll : iTuomil l201lF T As in Tuomi 
(2011a, 2011b), we also used the sample to assess 
the relative goodness of the different models by 
calculating their Bayesian probabilities given the 
data. For a positive detection of k+1 signals in 
the data, we required that the posterior probabil- 
ity of a k+1 -signal model was 150 times greater 
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than t hat of a k-signal m odel (e.g. iKass &: Rafterv 
199.4 iFeroz et al.ll2oTlh and that the amplitudes 
of all the signals were statistically distinguishable 
from zero. 

The Bayesian search we performed using differ- 
ent initial states in the parameter space as starting 
points of our posterior samplings, quickly located 
a signal with a period of 5.74 days. This signal 
satisfied our detection criteria since the probabil- 
ity of no signals is only 6xl0~ 8 , whereas a 1-planet 
model has a probability of essentially unity. The 
samplings did not reveal any other significant sig- 
nals in the data. Fig. [3] shows the probability den- 
sity function (PDF) for the returned periods as 
a histogram. The solid curve is a Gaussian pro- 
file with the same characteristics as the histogram, 
peaking at the measured period, and spread across 
a confined range of periods. The plot key high- 
lights some of the Gaussian parameters, in par- 
ticular the (T-width of the distribution shows how 
well constrained the signal is using this technique, 
as the range of periods here is only at the level of 
0.001 days. 

Fig. 0] shows a similar output PDF but for the 
possible range of eccentricities. This histogram 
does not have a Gaussian shape, more Poissonian, 
peaking around 0.10 or so. Again a Gaussian pro- 
file with parameters from the data is shown for 
reference. The parameters that describe the Gaus- 
sian are shown and the a is now at the level of 
0.07 in eccentricity. With a period of less than six 
days, putting the planet closer to the star than the 
tidal dissipation radius, we expect the actual curve 
is probably more circular, but with the relatively 
small number of data poi nts, we are seeing a more 
eccentric planetary orbit (jO'Toole et al.ll20081) . 

To test the plausibility of the non-zero ec- 
centricity model for t his system we ran the 
Lucv fc Sweenevi ( 197lh test on the data. The 
test returned a significance, p, of 6.7%, showing 
that the eccentric orbit is not significant at the 5% 
threshold level put forward by Lucy & Sweeney. 
Due to this result we also show a model fit that 
has the eccentricity fixed to zero in Table 01 

The final PDF we are interested in is shown 
in Fig. [5j which is for the semi-amplitude of 
the Doppler signal and shows, in accordance to 
our detection criterion, how it differs significantly 
from zero. The peak of the histogram is around 
6.0 ms _1 and the histogram again has a Gaussian 
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Fig. 3. — Histogram estimating the probability 
density function of periods from the Bayesian anal- 
ysis. Also overplotted by the solid curve is a Gaus- 
sian curve with the same mean and sigma as the 
distribution. 




-02i 3 : 0.69330502 



0.35634473 



0.4 



Fig. 4. — Histogram estimating the probabil- 
ity density function of eccentricities from the 
Bayesian analysis. Also overplotted by the solid 
curve is a Gaussian distribution with the same 
mean and sigma as the histogram. 
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shape, modeled by the solid curve and from the 
key we see that the a- width is less than 60cms _1 . 
From this semi-amplitude and estimated stellar 
mass we directly derive the minimum mass of the 
orbiting planet, finding a value of 15.9 M®. 

The final Keplerian curve describing this Doppler 
signal is shown in Fig. [6l phased to the 5.74 day 
period of the planet. It is clear that the data de- 
scribe the periodicity well. The uncertainties on 
the data points do not include estimated activity 
related jitter. Including the estimated jitter shown 
in Table [3] does not change any of the final orbital 
values since we treat the excess noise in the data 
as a nother free paramet er in our statistical model 



I see iTuomi et al.l 120111 ). Our sampling strategy 
of observing stars over consecutive nights helps 
to uncover these short periods, giving us efficient 
access to low-mass planets on short period orbits 
like HD773386. We also note that the 7 values we 
quote are with respect to the local mean of the 
data in Table. Q] 

One issue we have to contest with in this data 
set is the inter-dependency of the velocities we ob- 
served in a single night. On seven of the nights 
more than one radial velocity point was measured, 
meaning in the strictest sense, the data are not 
independent. However, to take this into consid- 
eration we have recently developed a model that 
employs an Auto- Regression and Moving Average 
(ARMA) algorithm to model the noise as an in- 
put prior into our full Bayesian model (for details 
and performance of the ARMA model we refer the 
reader to Tuomi et al. 2012, in preparation). We 
ran our ARMA Bayesian model on this data set 
and find the solution to be statistically similar to 
the ARMA free algorithm. We find that the mov- 
ing average components become statistically in- 
distinguishable from zero. Therefore, there are no 
correlations in the noise, at least not large enough 
for us to detect. 

5. Activity Indicators 

It has been well established that the CCF bisec- 
tor span (BIS) can correlate with radial velocity 
timeseries data, indicating that a given frequency 
detected in the velocities is not due to a genuine 
Doppler motion, but originates from the star itself, 
from processes that induce line asymmetries in the 
optical spectra of stars (e.g. iQueloz et al.l 12001 
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Fig. 5. — Histogram estimating the probability 
density function of semi-amplitudes for the data 
from the Bayesian analysis. Also overplotted by 
the solid curve is a Gaussian distribution with 
a matching mean and sigma to that of the his- 
togram. 
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Santos et al. 2010l) . Asymmetries can arise from 
spot patterns rotating across the disk of the star, 
blocking some of the light, and leading to mis- 
shapen spectral lines that correlate with the stellar 
rotation period. When fixed Gaussian profiles are 
fit to the final CCFs, they do not model the mis- 
shapen profile, and so the centroid of the Gaussian 
"wobbles" around the true value, in phase with the 
stars rotation period. 

One test of a true Dopplcr signal that is not 
related to stellar activity or velocity fields within 
the star, is to look for correlations between the BIS 
and the radial velocity timeseries. In the top panel 
of Fig. [7] we show the radial velocities against the 
BIS for HD77339. The solid line is the best linear 
fit to the data and reveals that there is no correla- 
tion between the two parameters. To place this on 
a statistical footing we also show the Pearson lin- 
ear correlation coefficient (r) within the plot, and 
the value of -0.08 shows that there is no statistical 
correlation between the radial velocities and the 
CCF BIS values, highlighting that the signal we 
detect is not due to line asymmetries in the star. 
Also a periodogram analysis does not reveal any 
power around a period of 5.7 days. 

As mentioned above, magnetically active re- 
gions (spots) rotating across the star disk due to 
the rotation of the star can lead to frequencies 
appearing in the radial velocities. Given we mea- 
sured the star to have a logi?HK of only -5.05, we 
would not expect the star to be very spotty and 
also we would not expect a rotational period as 
short as ^6 days. Also, our vsim measurement is 
very small (2.33±0.05 kms -1 ), confirming a longer 
rotational period is preferred for HD77338. In 
any case we still test if the chromospheric activ- 
ity indices correlate with the radial velocities for 
HD77338. 

The lower panel of Fig. [7] shows the HARPS 
activity S*-indices against the radial velocities. 
These indices were computed directly from the 
extracted ID HARPS-DRS spectra, following the 
methods we have previou s ly lai d out in detail in 
Jenkins et all (|2006l I2008L l201lh . The plot shows 
there is no apparent correlation between the two 
properties. The solid linear fit has only a small 
slope and we find the correlation coefficient, r, to 
have a value of -0.16, showing no significant rela- 
tionship between the activities and the radial ve- 
locities. Again, a periodogram analysis shows no 
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Fig. 7. — Correlations between the radial veloci- 
ties and the activity indicators that can be drawn 
from the HARPS reduced spectra. The top panel 
shows the CCF bisector span against the radial ve- 
locities and the lower panel shows the activity S 
indices for HD77338 as a function of the velocities. 
The solid lines show the best fit linear correlations 
between the parameters and the correlation coef- 
ficients are highlighted in the panels. 
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significant power anywhere across all periods. We 
do note that this fit is unweighted, and the lower 
activity point with low S/N and large uncertainties 
is probably pulling the slope lower than it should 
be given the rest of the data set. 

We note that it is unlikely that rotationally 
modulated spot patterns could give rise to the de- 
tected radial velocity variations since the signal 
has been stable for at least seven years. Also the 
measured rotational velocity is very low, mean- 
ing a rotation period of less than six days would 
require an inclination angle of ~18°, assuming a 
spherical rotating body. Such a low inclination 
has probability of only 4.9%. The very low chro- 
mospheric activity would also argue against a ro- 
tational period of less than six days, along with the 
presence of large spot patterns, and given that the 
star is of K spectral type, we expect a lower con- 
trast between the stellar photosphere and spots, in 
comparison to typical Sun-like G-dwarfs, meaning 
the radial velocity induced variation from any sta- 
ble spot patterns would be lower in comparison to 
more Sun-like stars. 

In addition to correlations between the BIS and 
S'-indices with the radial velocities, we also test for 
correlations between the CCF full-width at half 
maximum (FWHM) or the contrast of the CCF 
with the velocities. The FWHM and the contrast 
of the CCF has also been shown to be useful in- 
dicators of line variations, but at a lower level of 



signifi cance for these types of stars ([Boisse et al 
20111 ). Our correlation analyses found no strong 



correlations between these parameters and the ve- 
locities, strengthening the case for the reality of a 
Doppler signal. The final conclusion that can be 
drawn from analysis of the properties of the CCF 
and the chromospheric activity values is that the 
short period signal we have discovered in the radial 
velocity timeseries of HD77338 is from a genuine 
short period planet orbiting the star, with a low 
minimum mass that gives rise to the possibility 
of a planet with a non-negligible rocky mass frac- 
tion. The final parameters describing the orbit of 
HD773386 are shown in Table 0] for our best fit 
(top) and with a fixed eccentricity of zero (bot- 
tom). 



6. Transit Follow-up 

Short period planets orbiting bright stars can 
potentially yield transit dips in the light profile 
of their host stars, and by doing so, can lead to 
detailed follow-up of atmospheric physics, such as 



GJ12146 (e.g. 


Pont et al. 20081 Sw 


ain et al. 2009: 


Beaulieu et al. 


20101 iBean et al.l 


20111). Since 



HD77338 is bright (but not too bright to preclude 
transit detection using 1-m class telescopes e.g. 



Jehin etal 



(IS)), having a y-band magnitude 
of 8.63, and the planet has an orbital period of 
less than six days, giving rise to a transit probabil- 
ity of 6.4%, we performed photometric monitoring 
around the predicted transit time from the orbital 
solution to test if HD773386 could be added to this 
class of bright transiting planets. 

We performed the search for the transit of this 
planet using our current best fit solution for the 
system parameters, however we suspect a solution 
with a more circular orbit may be closer to the 
true nature of the HD77338 system. Therefore, 
assuming a circular orbit and an inclination angle 
of 90°, then we can estimate the transit time of 
the planet from 



P 

— arcsm 

7T 



(1) 



where P is the orbital period, R* and R p are 
the stellar and planetary radii respectively, and a 
again is the semimajor axis. Assuming a planetary 
radii commensurate with that of Uranus, then the 
transit time is predicted to be 2.9 hrs. Even if the 
planetary radius is upwards of a Jupiter radius, 
the transit duration is still estimated to be 2.9 hrs 
as it contributes such a small amount to the net 
radii in comparison to the star. 

The star was observed on the 25th of March, 
2012 and we used the lm-Cerro-Tololo Inter- 
American Observatory (CTIO) Telescope, along 
with the Y4KCam instrument, located at the 
CTIO observatory, in Directors Discretionary 
Time, to gain the necessary photometric precision 
to search for a transit of the discovered planet. 
Y4KCam is a 4064x4064 pixel CCD camera with 
a field-of-view of 20x20 arcmin 2 and a pixel scale 
of 0.289"/pixel. The standard readout time of the 
camera is 46s, which we reduce down to 16s by 
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Table 4: Maximum a posteriori orbital parameters and their 99% credibility intervals for HD773386 for 
unfixed (middle column) and fixed (right column) eccentricity. 

Parameter Model 1 Model 2 

Unfixed Eccentricity Fixed Eccentricity 
Orbital period P (days) 5.7361 [5.7346, 5.7376] 5.7361 [5.7345, 5.7374] 
Velocity amplitude K (m/s) 5.96 [4.22, 7.70] 5.96 [4.27, 7.62] 

Eccentricity e 0.09 [0, 0.34] 0.0 [fixed] 

lo (rad) 3.9 [0, 2tt] 0.0 [fixed] 

M (rad) 3.5 [0, 2tt] 5.8 [0, 2tt] 

m sin i(M e ) 15.9 [10.6, 20.6] 15.5 [11.1, 20.9] 

Semimajor axis a (AU) 0.0614 [0.0580, 0.0645] 0.0611 [0.0581, 0.0642] 
7 (m/s) -0.64 [-1.81, 0.53] -0.66 [-1.79, 0.36] 

rms (m/s) 1.74 1.71 

Nobs 32 32 

The uncertainties on the m sin i and semimajor axis consider the uncertainties on our stellar mass 
estimate. The 7 offset is the value after subtracting off the mean of the data set. 
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Fig. 8. — Photometric light curve of HD77338 covering the predicted time the planet would pass in front 
of the host star. The red curve is the model prediction for how long and how deep the transit should have 
been if it was observed. 



11 



employing a spatial binning of 2x2. 

We employed a strategy of observing the star 
at integration times of between 7-10 seconds each, 
with the telescope operating at heavy defocus to 
get PSFs with a typical radius of ^28 pixels. We 
obtained 786 data points with typical counts in 
the range 20'000-25'000, across a time baseline of 
over 6hrs, crossing the predicted center of tran- 
sit and the predicted egress of planet off the limb 
of HD77338. The observed images were then pro- 
cessed using a custom-made pipeline that performs 
the necessary trimming, bias subtraction, and flat- 
field correction to prepare the images for final 
analysis. 

Fig. |S] shows the normalised photometric light 
curve of HD77338 when the planet was estimated 
to pass in front of the star. We find no suggestion 
of a transit around the predicted time of center of 
transit of JD 2456012.527888. The rms scatter of 
our unbinned photometric data is 0.0028 meaning 
we could significantly detect any transit securely 
with a depth of only 0.008 percent below the con- 
tinuum level. 

For comparison, the red curve highlights the 
predicted depth and time of transit for HD773386 
given the orbital solution we found from the radial 
velocity data and a planetary radius commensu- 
rate with that of Saturn. Since no transit dip was 
detected we can say that either the inclination of 
the system is so low that the planet does not pass 
in front of HD77338 in our line of sight viewing 
angle, and/or the orbital solution we present is 
not secure enough to pin down a transit detec- 
tion. The first of these may be the case and is 
something we can not do anything about, however 
the second can be remedied with more radial ve- 
locity data, which we are still acquiring to search 
for additional planets in the system. The Lucy & 
Sweeney test we performed also indicated that a 
circular orbit can not be ruled out and therefore 
we should again search for a transit of HD773386 
using the circular orbital solution to test whether 
or not we can detect any occultation of the planet. 

7. Kinematic Motion 

The Hipparcos proper motion of this dwarf star 
is 41.13±0.67 mas in /i ra and -270.80±0.67 mas 
in /Zdoc, meaning it is classed as a high proper 
motion star. We used these values to mea- 



sure the three kinematic ve locity components 



(U,V,W) (I.Tenkins et al.1 l201lh 



and find veloci- 
-24.8±1.3 kms-\ 



ties of 39.1±2.0, -27.7±1.8, and 
respectively. 

Fig. [9] shows the distribution of all the stars in 
Jenkins et al. in the UV-plane, along with the 
position of HD77338 (filled red circle), the bound- 
aries of th e young (rec tangle), and old (ellipse) 
disks from lEggenl (jl969l ). HD77338 is found close 
to the boundary between the young and old disk. 
Fig.HOlshows a Toomre diagram for the same sam- 
ple and we see here that HD77338 is consistent 
with either the thin or thick disks, in agreement 
with the UV-plane kinematics. We also show the 
positions of all other hosts stars that have planets 
with minimum masses <0.05 Mj with the colours 
of the points representing the metallicity of the 
host star, from metal-poor being blue through to 
red at the metal-rich end. 

In comparison to the other low-mass planet 
hosts, HD77338 is fairly high in the energy plane 
and very central in the V kinematic velocity plane, 
along with being the most metal-rich. As ex- 
pected, most of the host stars are located around 
the region that contains mostly thin disk stars, 
with a few probable intermediate and/or thick 
disk hosts. However, the two stars HD4308 and 
HD20794 are found to have kinematics that place 
them well in the thick disk region and encroach- 
ing on the halo proxy perimeter. Also the blue 
colours of these two show they are metal-poor 
stars which fits in with their kinematics, whereas 
the hosts with kinematics that place them in the 
thin disk have a range of metallicities from really 
metal-poor through to super metal-rich. These 
two thick disk hosts have very eccentric orbits 
around th e galactic center (Appendix 1; code 
taken from Scholz et al. 1996) which helps to con- 
firm their membership of the thick disk since older 
stellar populations tend to have eccentric orbits 
and high velocities. For example, halo stars in 
general have very eccent ric orbits and attain ve- 
locities above 180 kms~ 1 ( Nissen fc Schuster 2010l ) 
in the Toomre diagram. If these stars are in- 
deed bonafide members of the thick disk then 
this helps to show that low-mass planets were 
probably forming early in the formation of the 
galaxy, if the thick disk was formed t hrough 
some hierarchal merger or heating process ( Statler 
119881: lAbadi et al]l2003t IVillalobos k HelmiH 2008: 
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Brook et al.l 120041 ) . However, other recent studies 
refute the e xistence of a thick disk being born in 



this fashion (jBovv et al.ll201lL 12012) 



To further test if HD77338 is also a thick disk 
member we perform an analysis of its orbital mo- 
tion through the galaxy. Appendix 2 shows the 
orbital plots for HD77338 (top four panels) and 
we see that it has a fairly eccentric orbit and it 
travels to ±0.2kpc above and below the plane of 
the galaxy. The lower panels show the same data 
for the Sun, a typical thin disk star, and we clearly 
see the circular orbit around the galaxy. Also, it 
only travels half of the distance in Z, above and be- 
low the galactic plane, in comparison to HD77338. 
These points may indicate that HD77338 is actu- 
ally a metal-rich member of the intermediate/thick 
disk, rather than the thin disk, but thin disk mem- 
bership can not be ruled out. 

The ratio of a elements in comparison to iron 
can allow one to distinguish between different stel- 
lar populations in some cases. To parameterize the 
a abundance of HD77338 we used the a elements 
shown in Table [5J In general the measured abun- 
dances depend on the adopted micrq turbulent ve- 



locity, V( fsee iPavlenko et alj 120121 ). Using the 



criteria of the absence of any dependence between 
log N(Fe i) and log N(Fe n) against the excitation 
potential E" to select the best absorption lines, we 
found the best fits to the observed profiles of the 
Fe I and Fe II ions with a V* = 0.5 km/s. This 
value was then adopted for the determination of 
all abundances listed in Table [2j 



[a/Fe] = [a/H] - [Fe/H] 



(2) 



We used the mean of the elements in Table [5] 
to obtain a logA(a) value of -5.74±0.09 dex for 
HD77338, and a solar value of -6.08±0.11 dex, giv- 
ing rise to a [a/H] of +0.34 dex. By application 
of Eq n [2]we get the final [a/Fe] value of -0.01 dex. 

We find no real enhancement of the a abun- 
dance for a star of such metallicity. iBensbv et al. 
([20051 ) show abundance trends of such elements 



and the values of our elements agree well with both 
the thin and thick disk stars for a super metal-rich 
star. In fact, Bensby et al. show there there is 
no significant enhancement of most a elements for 
thick disk metal-rich stars in comparison to metal- 
rich stars from the thin disk. The conclusion that 
can be drawn from the kinematics and abundance 
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Fig. 9. — A plot in UV space of the kinematic 
motion of all stars in Jenkins et al. (2011). The 
filled red point marks the position of HD77338. 
The solid box and solid ellipse mark the reg ions of 
the yo ung and old disk, respectively, from lEggenl 
(|l96flh . 
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Fig. 10. — A Toomre diagram for the sample of 
stars in Jenkins et al. (2011) (black points). The 
coloured star corresponds to HD77338 and the 
other coloured circles are all radial velocity de- 
tected planets with sub-Neptune masses and with 
colours ordered from the most metal-poor (blue) 
to the most metal-rich (red). The data was taken 
from the Exoplanet Database as of February 2012. 
The dashed lines mark regions of constant kinetic 
energy. 



13 



analyses is that HD77338 is most likely a thin disk 
star, or in an intermediate region between the thin 
and thick disks. 

7.1. Elements 

The [X/H] abundances for all the a elements 
we analysed are much higher than the solar values, 
in agreement with the super metal-rich nature of 
HD77338. The silicon abundance is significantly 
enhanced if we consider both the neutral and ion- 
ized silicon lines, however as can be seen from Ta- 
ble [2] the Si II analysis only employs two useful 
lines, in c omparison to 1 3 for th e Si I abundance. 
Recently, iTrevisan et al. (2011) also found a high 
silicon abundance for HD77338 of +0.47±0.05 dex, 
in good agreement with the value we find here. 

Given that silicon is enhanced in general for 
thin and thick disk stars in comparison to the 
amount of iron, it becomes a very important ingre- 
dient in planet formation models. Large amounts 
of silicon, and other elements, in the interiors 
of stars, and by association, the chemical make 
up of proto-planetary disks, mean a higher disk 
surface density in general. Changing metallicity, 
and hence disk surface density, can have a pro- 
found impact on planet formation affecting both 
the masses of planets that can f orm and also by 



how much they can migrate (see iMordasini et al 



20121 ). Metal- rich systems like HD77338 generally 
give rise to gas giant planets, not lower-mass plan- 
ets like HD773386. However, population synthesis 
modeling is a statistical approach and gives gen- 
eral outcomes for the observed population. It docs 
not preclude the formation of very low-mass plan- 
ets in metal-rich systems, but rather concludes 
that higher mass planets are more likely around 
these stars in comparison to stars with a lower 
iron abundance. 

The other elements are all found to have values 
of [X/Fe] commensurate with that of the iron in 
HD77338. Indeed, iron has long been used as a 
strong proxy for the overall metallicity ([M/H]) of 
stars like the Sun. Therefore, for this system we 
might expect that any core in H D773386 is e l ement 
hea vy and possibly silica te rich. ISato et al.l (|2005i ) 
and lFortnev et al] ([20061) show that the transiting 
planet HD 1490266 must be around 2/3 heavy el- 
ements by mass, and given that the star is super 
metal-rich, this points to a connection between the 
abundances of elements in stellar atmospheres and 



those of their retinue of planets. A further system- 
atic search for the transit of HD773386 could help 
to shed light on this issue. 

8. Mass and Metallicity Discussion 

HD77338 is one of the most metal-rich stars 
currently k nown to host a plane t. If we take 
the recent Trevisan et al. ( 201 lh [Fe/H] value 
of +0.47±0.05 dex for HD77338 then only one 
planet host star has a metallicity higher than this, 
HD126614 A with a metallicity of +0.56±0.04 dex 
([Howard et al.ll2010h . However, HD126614 A has 



an M dwarf companion at only 33 AU separation 
from the host star, meaning HD773386 would be 
the most metal-rich single star known to host an 
exoplanet. Also, the low-mass nature of HD773386 
(m sin i < 0.05Mj) helps to populate the low-mass 
metal-rich bin and shows that metal-rich, very 
low-mass planets may be plentiful. 

The high abundance of all the elements we 
have analysed also indicates that the proto- 
planetary disk left over from the formation of 
HD77338 was rich in planet building material. 
The established preponderance fo r gas giant plan- 
ets to favour metal-ric h stars (jGonzalez 1997 ; 
Fischer fc Valentil l2005t ISousa et all l201ll) can 
be seen in Fig. QTJ We also see that metal-rich 
stars tend to cover the entire phase space of plan- 
etary masses. Above a host star metallicity of 
~0.2 dex there are planets covering the whole 
regime from low-masses to high-masses, in dense 
clusters. However, for lower metallicity systems, 
particularly at sub-solar metallicities, there are 
regions free from any planet detections, or at least 
less densely packed. This shows why metal-rich 
stars are so highly prized in the hunt for exoplan- 
ets and for better understanding the nature of 
planet formation and migration. 

8.1. Low-mass Systems 

Lower mass systems may show a somewhat dif- 
ferent met allicity distribution in comparison t o 



gas giants ([Udrv et al.l 120061 iNeves et al.l l2009h . 



with no clear metallicity bias. Models indicate 
that gas giant planets form outside the ice line 
boundary, taken to be ^5 AU, in metal-poor disks 
( Mordasini et al 1 l2012h meaning smaller rocky/icy 
planets can form interior to this boundary, giv- 
ing rise to the observed population of short period 
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Table 5: Orbital and stellar parameters for all radial velocity detected host stars with sub-Neptune mass 
planets. 

Planet m sin i Somimajor Axis Period eccentricity K M* [Fc/H] V SpT U,V,W 



M j AU days ms AI q dex mags kms 
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Fig. 11. — Plot of radial velocity detected plan- 
ets host star metallicity against minimum mass 
of the planet, taken from the Exoplanet Database 
(jWright et al.l 120111 ). Only planets around dwarf 
stars were included. 
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Fig. 12. — Metallicity against minimum mass 
for all radial velocity detected planets with 
a minimum mass less than that of Neptune 
(m sin i < 0.05 Mj). The horizontal dotted line 
highlights the canonical boundary where runaway 
gas accretion onto the growing core is expected to 
occur. The dashed line marks the lower bound- 
ary in mass and the data points have been scaled 
in size by their orbital period. The filled circles 
connected by a straight line show the position of 
HD77 3386 using our metall icity and the metallic- 
ity of Trevisan et al.l ( 201ll ). respectively. 



rocky planets around stars with a lower metal con- 
tent. 

Fig. H2] shows the distribution of sub-Neptune 
mass planets (m sin z<0.05Mj) as a function of 
metallicity using data taken from the Exoplanet 
Database! as of February 2012. All values we 
have used are listed in Table El One of the most 
striking features we see is the lack of planets with 
high metallicities and very low masses, in the lower 
right corner of the plot. Even though the num- 
bers are relatively small at present, there appears 
a lower boundary that increases with metallicity, 
which would be contrary to the presumed biases in 
current radial velocity surveys given the number of 
metal-rich programs currently running. HD77338 
also appears to follow this trend given that it has 
a very high metallicity and a minimum mass com- 
mensurate with that of Uranus in our solar system. 

A strong bias that could be manifest here is 
from the distribution of the orbital periods of the 
sample. Given the high fraction of short period 
gas giants in metal-rich systems, it is possible that 
lower mass planets reside at longer orbital periods, 
meaning they will be biased against due to the well 
known radial velocity bias agai nst longer perio d 
and low amplitude objects (see Cumming 2004) . 
To test this we scale the data points in Fig. [12] by 
the ratio between each planet's period and that 
of the longest period planet in the set, GJ876e. 
This shows us that this is not the case and that 
as expected, the planets are mostly short period 
planets at the perimeter of this boundary. It is 
also noticeable that beyond the cluster of plan- 
ets near the perimeter of the boundary we mark 
by the dashed line, there is a possible valley in 
the metallicity plane, then the rest of the sam- 
ple is found. A metallicity-mass correlation like 
this, at low-masses, would fit well with the core 
accretion model for planet formation, where the 
low-mass planets around super metal-rich stars 
attain a higher mass due to the abundance of 



planet forming material in t he disk (jlda fc Lin 



2004ri iMordasini et al.M2012l ). What is not ob- 
vious is why there would be a lower boundary. 
Possibly the planetesimals quickly accrete mate- 
rial and grow much faster than at lower metallici- 
ties, meaning the boundary traces the planet for- 
mation timescales for a given disk metallicity. In 



2 http: / /www. exoplanets.org 
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any case, if there is a lower boundary that holds 
to long period orbits for super metal-rich stars, 
this has consequences for the fraction of habitable 
Earth-mass planets in the galaxy, as super metal- 
rich stars will not be abundant in such planets. 

8.1.1. Monte Carlo Test 

To test if the lack of planets in the lower right 
corner of the metallicity-mass plane is statistically 
significant, and hence a metallicity-mass correla- 
tion for low-mass rocky planets is statistically sig- 
nificant, we perform a Monte Carlo analysis. We 
setup the model in two parts where we first gener- 
ate a random realisation of the sample of 36 low- 
mass planets we have tested in Fig. [T2j To gen- 
erate the masses in a robust fashion we take the 
observationall y driven mass distr i bution for lower- 
mass planets ( Butler et al.l 1200(1 ICumming et al.l 
2008). Lopez & Jenkins (2012, in preparation) 
show that a power law model such as this explains 
the observed turnover of the mass distribution at 
the lowest masses and hence the model is a robust 
representation of the current distribution of exo- 
planets. Eq n [3]shows the form of the mass function 
we consider. 



dN/dM = kM 



-1.2 



(3) 



Here dN/dM represents the frequency of plan- 
ets as a function of mass, M is the mass bin consid- 
ered, and k is a scaling constant to fit the observed 
population. We can then integrate this mass func- 
tion equation (Eq 11 |4|) to get the result shown in 
Eq n [5l and by taking the inverse we can generate 
the cumulative distribution function (CDF; Eq n [S]) 
that can be used to draw masses randomly from 
the observed population of exoplanets. 



dN = k 



A I 



N 



M- 12 dM 



-5k 



M ' 



(4) 



(5) 



Now if we set the constant k equal to -0.014 
this will normalize the function such that we have 
a probability density function and the correspond- 
ing CDF is given by: 



CDF(M) 



0.069 

N 



The CDF allows us to draw random masses for 
planets and restrict the values to be within the 
range between 0.006 to 0.050 Mj to match our 
region of interest. We can then randomly gener- 
ate metallicities, and in this test we simply draw 
from a uniform distribution of [Fe/H] values in the 
range from -0.5 to +0.5 dex, which covers all the 
values of metallicities in our test. 

Once we have built a random sample of plan- 
ets we can then add the boundary and test how 
many times there are planets that reside under 
the boundary region we identify in Fig. [121 We 
run the code l'000'OOO times to ensure a high 
level of robustness in the final probability mea- 
surement. Our test reveals that for a sample of 
36 stars, assuming a uniform metallicity distribu- 
tion, 99.9993% of the time there are planets to be 
found below the boundary region. Therefore, this 
low-mass planet desert is statistically significant 
at almost the ~4.5<r level, under these test con- 
ditions. For a sample of only 10 planets we still 
find that 96% of the time there are planets below 
the boundary region, and for a sample of 5 planets 
we find a percentage of 81%. These tests indicate 
that there is some significant correlation between 
metallicity and mass for low-mass planets, at least 
in the metal-rich regime, that would be important 
to quantify in the future. 

We do note that our results can vary due to the 
assumed metallicity model that we use to draw the 
sample of random metallicities from. For instance, 
the result will become more significant for a distri- 
bution that follows the current observed distribu- 



tion for more massive exoplanets (see ISousa et al 



2011), and vice- versa for a distribution shaped in 



the opposite fashion. 

As for the possible valley, the sample is as yet 
too small to draw statistically significant conclu- 
sions and therefore difficult to test without de- 
tailed modeling of the core accretion method of 
forming planets, therefore we are only pointing out 
the possibility that there are two classes of low- 
mass planets that have a metallicity dependency. 
In particular, magnitude or distance limited sam- 
ples will be biased towards sub-solar metallici- 
ties as they will be governed by the metallic- 
ity dis tribution of the local galactic neighbour- 
hood ( Holmberg et al. 20091) . However, to fully 



(6) 



probe these samples it is necessary to simulate 
the radial velocity data for each star individu- 
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ally (|0 'Toole et alj2009[ ) to better understand the 
completeness of the bins we have discussed. 

8.1.2. What does Kepler Say? 

The Kepler Space Telescope has recently re- 
leased a plethora of transiting planet candidates, 
that range in radius from bel ow Earth radii up int o 
the gas giant planet arena ([Borucki et al.l 1201 lh . 
Follow-up of Kepler targets are difficult due to 
the relative faintness of the candidate stars, how- 



ever, iBuchhave et al.l (J2012J) have published the 
first large set of spectroscopic metallicity for 152 
Kepler planet-host stars that can be used to draw 
statistically significant conclusions on Kepler host 
star abundance patterns. 

In Fig. Q2] we show the Kepler results published 
in Buchhave et al. in the metallicity - radius plane. 
We only show the stars that host planets with radii 
of 4i? e or less, and we do not distinguish between 
multiplanet systems and single planet host stars. 
The 4i?0 limit was chosen simply because typi- 



cal mass-to-radius re lationships, e.g. M p = R p 



2.06 



(jLissauer et al.l 120111 ). would give rise to around 
the Neptune-mass limit highlighted in Fig. [T2] 
The paucity of planets we find in the metallicity 
- mass plane from the radial velocity data could 
therefore be manifest in the Kepler data too. How- 
ever, we note that a single mass-to-radius relation- 
ship is surely unrealistic due to the possible diver- 
sity of rocky planets (see lSeager et al.ll2007l ). 

Studying the lower right corner of Fig. [T3J the 
high metallicity and low mass regime, we do find a 
paucity of planets when we compared to the lower 
left corner, or the low metallicity and low mass re- 
gion. To better highlight this relative planetdesert 
we bound the region by a dashed box and also a 
similar rising boundary line (solid curve) as we 
show in the metallicity - mass plane. The Kepler 
data is in good agreement with the radial velocity 
sample in respect to this paucity of planets, al- 
though not likely in absolute value, with only two 
planets significantly below the boundary region we 
highlight. In fact, within the boxed region, pa- 
rameterised by the host star having a metallicity 
of solar or above and the transiting planet having 
an Earth radius or lower, there are only 8 plan- 
ets, and only one of these planets is currently in 
a single system. Therefore, it could be that really 
low-mass planets can exist in the super metal-rich 
regime if they are part of multi-planet systems, 



since under core accretion the planetesimals will 
compete for material to form and this will mean 
less material for each of the planetesimals to reach 
higher core masses. 

If these features however are real, then for a 
given stellar/disk metallicity it appears that plan- 
etesimals quickly grow to their boundary region 
and stop there, or they then quickly make the tran- 
sition to around 2 to 3 times their boundary mass 
and either continue growing to higher mass plan- 
ets before the disk dissipates, or stay where they 
are due to depletion of the disk material. 
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Fig. 13. — Kepler results for 152 planet-hosting 
stars detected by Kepler in the metallicity - radius 
plane. The dashed box region in the lower right 
corner highlights the lack of low-mass Kepler tran- 
sits in comparison to the more metal-poor region. 
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A. Orbital Plots for HD4308 and HD20794 

B. Orbital Plots for HD77338 and the Sun 
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Fig. 15. — The calculated orbital motion of HD77338 through the galaxy in all three spatial dimensions X, 
Y, Z, along with the distance from the galactic center are show as a function of time are shown in the top 
four panels, from top left to lower right, respectively. The same orbital calculations are shown in the lower 
four panels for the orbit of the Sun around the galaxy. 
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